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Abstract 

We employ the point canonical transformation (PCT) to solve the D-dimensional Schrodinger 
equation with position-dependent effective mass (PDEM) function for two molecular pseudohar- 
monic and modified Kratzer (Mie-type) potentials. In mapping the transformed exactly solvable 
D-dimensional [D > 2) Schrodinger equation with constant mass into the effective mass equa- 
tion by employing a proper transformation, the exact bound state solutions including the energy 
eigenvalues and corresponding wave functions are derived. The well-known pseudoharmonic and 
modified Kratzer exact eigenstates of various dimensionality is manifested. 
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I. INTRODUCTION 



The solution of the Schrodinger equation with position-dependent effective mass (PDEM) 
for an arbitrary central potential has attracted attention over the past years (cf., e.g. [1] and 
the references therein) . The motivation in this direction arises from considerable applications 
in the different fields of the material science and condensed matter physics. For instance, such 
applications in the case of the bound states in quantum system [2], the nonrelativistic Green's 
function for quantum systems with the position-dependent mass [3] , the Dirac equation with 
position-dependent mass in the Coulomb field [3], electronic properties of semiconductors 
[4] , ^He cluster [5] , quantum dots [6] , quantum liquids [7] , graded alloys and semiconductor 
heterostructures [8,9], the dependence of energy gap on magnetic field in semiconductor 
nano-scale quantum rings [10], the solid state problems with the Dirac equation [11], etc. 
Almost all of those works mentioned above were focused on obtaining the energy eigenvalues 
and the potential function for the given quantum system with the PDEM function. The 
wave functions were either obtained by the solutions to the Schrodinger equation with the 
constant mass, or a few lower excited states were obtained by acting of the creation operator 
on the ground state. The effective potentials are the sum of the real potential form and 
the modification terms emerged from the location dependence of the effective mass [12]. 
Taking into consideration the PDEM, the main concern, in this work, is in obtaining the 
energy spectra and/or wavefunctions of the D-dimensional Schrodinger equation with a given 
PDEM for central potentials by the point canonical transformation (PCT) approach [1] . In 
modern theory, high dimensions are also of interests in many fields [13-25]. 

In this work, we employ the PCT to solve the D-dimensional Schrodinger equation with 
PDEM for the pseudoharmonic [26-31] and modified Kratzer [32-37] potentials through 
mapping this wave equation into the well-known exactly solvable D-dimensional Schrodinger 
equation with constant mass for a given PDEM function [1]. Indeed, the PCT approach 
has enabled us to obtain the exact effective mass bound state solutions including the energy 
spectrum and corresponding wave functions in any dimension for the exactly solvable classes 
of quantum molecular potentials. 

This work is organized as follows: in section 2, we introduce the methodology. Section 
3 is mainly devoted to obtain the exact bound state energy eigenvalue and eigen function 
solutions of the D-dimensional Schrodinger equation with a given PDEM function for two 
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exactly solvable diatomic molecular potentials. We give a few concluding remarks in section 
4. 



II. METHODOLOGY 

The D-dimensional PDEM Schrodinger equation for a central potential V{r) takes the 
form 

where the position-dependent mass distribution m(r) is a real function. Here the wave 
functions (x.) belong to the energy eigenvalues E and V{r) stands for the D- 

dimensional standard central potential in the configuration space coordinates. Here r repre- 
sents the D-dimensional radius ( Yli=i ^1 ) ■ Atomic units will be used throughout, with 
h/27r — h — 777-0 — e — 1. Going over to a spherical coordinate system with D — 1 angular 
variables and one radial coordinate we can write [21-25] 

where y^^'^.^^ ,^ (x) represents contribution from the hyperspherical harmonics that arise in 
higher dimensions with x represents the D-dimensional position vector. With the substitu- 
tions 



D 



^v„^<;— ;'(x)) = (v„^) . (v„<— + ^v?,w;-:?(x). 



^2 _ {D-1) d l{l + D- 2) 



= \ 

1 d f . D-2n 9 \ Ll_, 

sm Bd^i- 



(3) 



.sin^-^^D-i-^^i^-i V dOD-iJ sin^^i,.! 

and the wave functions (2) into Eq. (1) will result in the following time independent D 

dimensional PDEM Zpth partial wave radial Schrodinger wave equation 
£ m'(r)(D -II d\ Id(Id + 1) 



+ 



+ 2m{r)[E-V{r)]^Rr,,i{r)^0, (4) 



dr^ m{r) \ 2 r dr J 

where m'{r) = dm{r)/dr and the transformation I = I + {D — 3)/2, D > 2 and V|, 

stand for the D-dimensional angular momentum and Laplacian, respectively. Moreover, the 
D — 1 can be obtained through inserting I — 1 and or Id — and —1. In the present work. 
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we are concerned in bound states, i.e., E <0. We also should be careful about the behavior 
of the wave function R{r) near r — and r — > oo. It may be mentioned that Ri{r) should 
be normalizable [8]. Obviously, it can be shown, most readily, that for a constant mass, i.e., 
m'{r) = case, the above equation with potential function V{s), angular momentum A(/) 
and energy spectrum E reduces to the usual case of independent-position mass: [19,21-23]: 

where Ad(/) = A(/) + {D — 3)/2. Consequently, the solutions for a particular central po- 
tential V{s) are the same as long as D + 2A{1) remains unaltered. For example, the s-wave 
eigensolution (Rq) and energy spectrum (E) in four-dimensional solutions are identical to 
the p-wave two-dimensional solutions (A(Z) — O^D — 4) — *-(A(Z) — 1, D — 2). For more de- 
tail on inter-dimensional degeneracies the reader may refer to, e.g. [19,23,38,39]. We apply 
the following point canonical transformation (PCT) s — >■ r (i.e., mapping function s = q{r)) 
with the substitution of the wave function 

Rn,i[r) = — , (6) 

into Eq. (5) and after some simple algebra, the transformed Schrodinger equation becomes 

dr'^ \ g q' ) dr \g q' g j 

- Ad{1) (A^(0 + 1) (^^) + 2 {qf [e - V{q{r))\ | i?„,A(o(r) = 0, (7) 

where the primes denote differentiation with respect to r. Further, comparing Eq. (7) with 
Eq. (4), we find the following PCT transformations: 

1 «' = m(,V(r) (8) 



2m(r) \ q' g m{r) 



V{r) = ^,V{s), (9) 
m[r) 

En = ^En, (10) 

m[r) 



and 



AMD (A.(0 + 1) (9 ^ '-^^ - - I m^nim - FW)l . (11, 

where F{x) — ^ — \ (^) ■ Therefore, Eq. (6) and Eqs. (8)-(ll) can be also found for any 
potential system with PDEM. 
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III. APPLICATIONS 



We solve the D-dimensional PDEM Schrodinger equation exactly for two potentials: the 
pseudoharmonic potential [26-31] and the modified Kratzer molecular potential [32-37] . The 
transformation function g{r) will be found for the given effective mass function m(r) = rrior^ 
and the selected PCT function q{r) — . 



A. pseudoharmonic potential 

The pseudoharmonic potential is given by [26-31] 

as a reference potential where V"e = \ht1 is the dissociation energy between two atoms in 
a solid with k is the force constant and Tg is the equilibrium bond length. The eigenvalue 
problem (5) for V{s) in (12) can be solved analytically to get the exact D-dimensional results 
for energy eigenvalues and wavefunctions of this system (in units in which mo — h — 1) as 
[26-30]: 

KmD = -2^e + (4n + 2 + V(^ + 2A(0-2)2 + 8Kr2) , (13) 



and 



Rn,A{l) 



xF{-n, V(A(0 + D/2 - 1)^ + 2Kr2 + 1-J^s'], (14) 



with 



2 (v^//-,)^^^"^"''''"'^'+'''''+' //! 
~ \ r{n + v/(A(iy + D/2 - 1)2 + 2K^ + 1) 
where n — 0, 1, 2, • • • and / — 0, 1, 2, • • • signify the usual radial and angular momentum 
quantum numbers, respectively and An^i being the normalization constant. Here F{—n, k + 
l,x) is a confluent hypcrgcomctric function. We follow Rcf. [1] by taking the power law 
PDEM function m(r) = mor^ and the PCT function q{r) = r", where mo is the rest mass 
and A and u are two non-zero real parameters. We consider only the case where {q')"^ /m{r) 
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is constant for which i/ = l + A/2, A7^2to avoid position-dependent energy. Inserting them 
into Eq. (6) and Eqs. (9)-(ll), we obtain 



^0 / ^l+A/2 



-V 



^l+A/2 I 



En,l — 



2 + A r 



-rjrl 



+ l + 2n+ v^(A(iyTD72^^1)2T^] C, 



(15) 
(16) 



and 



R^,{r) - a, (;^r)(^+^)(^/(^«+^/^-^)^+^^'^^+^)+^ exp (-|r^+^) 
X F (^-n, V(A(/) + D/2 - 1)^ + r]hi + 1; r/r2+^) , 



(17) 



where 77 = = 2moC/(2 + A) with C — vq/rriQ is a real potential parameter and A(Z) is 

(18) 



A(0 



{D-2)^ 1 



(L> + 2/ - 2)^ + (2 + A)^ - 2 (2 + AD). 



2 ' (2 + A) 

The radial wave function (17) must vanish as r ^ and r ^ 00. For the trivial case where 
A = 0, i.e., constant mass (m = mo), we obtain A(/) = / from Eq. (18). Hence, we find out 
that Eqs. (15)-(17) agree with Eqs. (12)-(14). The same procedure leads to the other exact 
solvable classes belong to various values of parameter A. For example, if we insert A = 2, the 
following potential function with its energy spectrum and wave functions are obtained 



V{r) 



2rf 

mo 



2 

r 1 



2\ 2 



n,l 



4:1] 

mo 



-rjrl + 1 + 2n + V(Ai(0 + £'/2 - 1)^ + TyVf] , 



(19) 
(20) 



and 



R.Jr) 



5^ (^,)V(2Ai(0+i^-2)^+4,V4+| |^_|^4j 



X F (^-n, V(A(0 + D/2 - 1)2 + jyVf + 1; rjr^'^ , 



with 



Ml) = -^-^j^ + y{D + 2l-2f-A{D-3), 



(21) 



(22) 



where n,l = 0, 1, 2, • • • . The pseudoharmonic potential can be treated exactly in three as 
well as in one and two dimensions. We note the special cases D = 1,2 and 3. For D = 2, 
with the customary notation / = M and r = p : 

,2 



En. 



M 



mo 



2 4M2 + A2 
-rjr^ + 1 + 2n + ^ / —3- + r| 



(2 + A)^ 



(23) 
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and inserting A = 



Er, 



mo 



1 



-rjrl + l + 2n+ -^l^^rfr^ 



and 



(33) 



(34) 



where n = 0, 1, 2, ■ ■ • . 



B. modified Kratzer molecular potential 

The modified Kratzer potential is [32-37] 

n.) = V;(^y, (35) 

as a reference potential. The eigenvalue problem (5) for V{s) in (35) can be solved analyti- 
cally to get the exact D-dimensional results for energy eigenvalues and wavefunctions of this 
system (in units in which mo — h—\) d& [32-37] : 

^n,A(0 = K ^ n = i = 0, 1, 2, • • • (36) 



2a2 (^1 + 2n + ^[D + 2k{l) -2f ^W^rl^ 



and 



X F (^-n, 1 + -\/[D + 2A(/) - if + S^r^; 2A;sj , (37) 

with a — 1/ (4T4re) and is the normalization factor. The wave number k for the modified 
Kratzer (pseudo-Coulomb) spectrum under consideration is 

A; = ^ (38) 



The PDEM function m(r) = m^r^ and the PCT function q{r) = are same as before. 
Inserting them into Eq. (6) and Eqs. (9)-(ll), we obtain 

/ l+A/2 _ ^ \ 2 

nr)=P^:T^ , (39) 
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EnJ = P- 



1 + 2n) (2 + A) + \ {2 + A)' [D + 2A(/) - 2]' + 32mor2p 



2 ' 



(40) 



and 



X 



(41) 



where P = 

is a real potential parameter and 
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IGmorgP 



(2 + A) 



(1 + 2n) (2 + A) + J (2 + A)' [Z^ + 2A(/) - 2]' + 32mor2p 



which are identical with those given in Refs. [32,34-36] when A is set to zero. The new angular 
momentum A(Z) is as given in Eq. (18). Particularly, setting A = with P — Ve/mo into 
Eqs. (39)-(41),we recover the constant mass Eqs. (35)-(37). This manifests the generality of 
our solution for the mass function given by m(r) = tuqt^. We note the special cases D = 1, 2 
and 3. For D — 2, with customary notation I — M, we obtain 



E., 



32ni.r-P- 



n,M 



[1 + 2n) (2 + A) + ^/T6WT32^n^ffPTV^ 



2 ' 



(42) 



and 



Rn,M[r) ^ hir^^ V e ; 4exp(-7ir^2l 

F {-n, 1 + + SVerl + 4A2; 27ir^+^) , 



with 



7i 



(2 + A)^ 1 + 2n + yi6M2 + 8yer2 + 4A2' 
wher n, M = 0, 1, 2, • • • . Further, when A = 

Smor^P^ 



En,M — P 



1 + 2n + ^4M2 + SmorlP 
Rn,Mir) = biA+V^^^^+^e exp (-7ir) F (^-n, 1 + 2^471^2^21^; 27ir) , 



where P = Vg/mo is a real potential parameter and 

AmoTeP 



7i 



1 + 2n + yi6M2 + 814^2' 



(43) 
(44) 

(45) 
(46) 

(47) 



which are consistent with those given in Ref . [28] . For D — 3 : 

32mor2p2 



EnJ - P - 



(1 + 2n) (2 + A) + 2^(1 + 2lf + SmorjP + A(A - 2) 



2 ' 



(48) 



and 



Rn,l{^) =Oir^ 'V (2+A) 



cxp [ — 72r^'''2 ^ 



72 



(2 + A)^ 
IQmoreP 1 



(49) 



2 + A)' 1 + 2n + 2,/(^I±l2§pr7^' 



where n, Z = 0, 1, 2, • • • . Interestingly for A = the above results in (48)-(49) reproduce the 
well-known three-dimensional relations obtained recently by Refs. [32,34-36,39]. The energy 
eigenvalues and wave functions are given in the form [32,34-36,40] 

8mor2p2 



En,l = P 



1 + 2n + ^(1 + 2lf + SmorlP^ 



(50) 



Rn,i{r) = &zr^v/8Fer-?+(2i+i)2 p {21 + if + Wevl; 2^2r^ , (51) 

(52) 



72 = 



AmoTeP 



l + 2n+ ^J{2l + lf + 8Verl 



For D — 1 (s-wave) : 



En^P 



32mor2p2 



;i + 2n) (2 + A) + 2^(1 + A)2 + Smor^P 



2 ' 



(53) 



and 



16moreP 



'(1 + A)^ 



(2 + A)^ 



(54) 



73 = 



where n = 0, 1, 2, • • • . When A = : 

8morgP2 



P. = P - 



(l + 2n + + 8mor2p) 



2 ' 



(55) 
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and 

Rn{r) = ft^r^Vs^^exp {-^sr) F (-n, 1 + 71+81^; 2^sr^ , (56) 

4moreP 
1 + 2n + yiTSK^ 

IV. COCLUDING REMARKS 

In this work, we have studied the exact PDEM Schrodinger equation in D-dimension for 
two diatomic molecular potentials, namely, pscudoharmonic potential and modified Kratzer 
potential. The exactly solvable constant mass D-dimensional Schrodinger equation has been 
transformed into the form similar to the effective mass by means of a proper PCT function. 
The mapping of the resulting transformed equation with the original effective mass equation 
provide us the required energy spectrum and wave functions for the potential system under 
consideration. We have applied this methodology to obtain the energy eigenvalues and the 
corresponding eigenfunctions of the modified Kratzer and the pscudoharmonic potentials. 
The exact bound state solutions of the constant mass Schrodinger equation for the pscu- 
doharmonic and modified Kratzer problems are recovered for various dimensionality upon 
inserting D — 1,2 and 3 and A = 0. 
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